Two-Dimensional Turbulence in the Inverse Cascade 
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q ■ Abstract. Numerical and physical experiments on the forced two-dimensional Navier- 

Q\ . Stokes equations show that transverse velocity differences are described by "normal" Kol- 



mogorov scaling < (Av ) 2n >oc r 3 and obey a gaussian statistics. Since the non-trivial scal- 
ing is a sign of a strong non-linearity of the problem, these two results seem to contradict each 



O. 

other. A theory, explaining this result is presented in this paper. Strong time-dependence of 



the large-scale features of the flow (u 2 oc t) results in decoupling of the large-scale dynam- 
ics from statistically steady-state small- scale random processes. This time-dependence is 



O 
> 

also a reason for the localness of the pressure-gradient terms in the equations governing the 
small-scale velocity difference PDF's. The derived self-consistent expression for the pressure 
gradient contributions lead to the conclusion that the small-scale transverse velocity differ- 
ences are governed by a linear Langevin-like equation, strirred by a non-local, universal, 
solution-depending gaussian random force. This explains the experimentally observed gaus- 
sian statistics of transverse velocity differences and their Kolmogorov scaling.. The solution 
for the PDF of longitudinal velocity differences is based on a smallness of the energy flux in 
two-dimensional turbulence. The theory makes a few quantitative predictions which can be 
tested experimentally. 
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1 Introduction 



Theoretical prediction of two inertial ranges, consequence of both energy and enstrophy 
concervation laws by the two- dimensional Euler equations, was and still is one of the most 
remarkable achievements of statistical hydrodynamics [1]. A direct and most important 
outcome of these conservation laws is the fact that if a fluid is stirred by a random (or 
non-random) forcing, acting on a scale If = 1/kf, the produced energy is spent on creation 
of the large-scale (/>//) flow which cannot be dissipated in the limit of the large Reynlds 
number v — > 0. This means that the dissipation terms are irrelevant in the inverse cascade 
range. Since the dissipation contributions are one of the most difficult obstacles on the road 
toward turbulence theory (see below), one can hope that in two dimensions the situation 
is greatly simplified. This hope is supported by recent numerical and physical experiments 
showing that as long as the integral scale Lj oc t% is much smaller than the size of the 
system, the velocity field at the scales Lj >> I >> If is a stationary close-to-gaussian 
process characterized by the structure functions 



where the pumping rate P is defined below [2] -[4]. Moreover, both numerical and physical 
experiments were not accurate enough to measure 



which were too small. This means that the observed probability density P(Au) was very 
close to symmetric one. This experimental fact differs from the outcome of the measurements 
in three dimensions where s n 's are very large when n is not small. Thus, the absence of strong 
(may be any) intermittency in two-dimensional turbulence and proximity of the statistics of 
velocity field to gaussian makes the problem seem tractable. 
The equations of motion are (density p = 1): 



S n = (u(x + r) - u(x)) n = (Am)™ oc (Pr) * 



(1) 



S2n+1 




« 1 



(2) 



d t Vi + VjdjVi = -dip + uV 2 Vi + fi 



(3) 
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and 

diVi = (4) 

where f is a forcing function mimicking the large-scale turbulence production mechanism 
and in a statistically steady state the mean pumping rate P = f ■ v. In the inverse cascade 
range the dissipation terms in (3) will be irrelevant. Neglecting it and multiplying (3 ) by V; L 
we obtain readily 

E = -v* = Pt 
2 

Thus, in this case the energy linearly grows with time. 
In this paper we define the force correlation function as: 

< /,(k)/,(k') >oc P{5 1] - |)^(k + k')^-i') (6) 
so that 



(5) 



(f(x + r)- fir-)) 2 oc P(l - Cos(k f r)) (7) 

It will be clear below that the forcing term enters the equations for the probability density 
of velocity differences exclusively through the expression (7) and in the limit fc/r << 1 it 
contribution is 0((/c/r) 2 ) which is a well-known fact. In the energy range we are interested in 
this work k/r » 1 and the oscillating contribution can be neglected leading to disappearence 
of the forcing scale from equation for the PDF. Thus the general expression for the structure 
functions is: 

S n (r) oc (Pr)f(^) f » (8) 
where the exponents 5 n denote possible deviations from the Kolmogorov scaling. If a statisti- 

3 

cally steady state exist in the limit L t » I » If, then all S n = since Lj oc £2. This would 
be prooof of " normal" (Kolmogorov) scaling in the inverse cascade range, provided one can 
show that the PDF P(Au) in the inertial range is independent on its counterpart in the 
interval I ~ If. This is the subject of the present paper which is organized as follows. In the 
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next Section the equations for the generating functions are introduced. Section 3 is devoted 
to a short analysis of the Polyakov theory of Burgers turbulence some aspects of which are 
used in this paper. Some physical considerations, which are basis for the developing theory, 
are presented in Section 4. S In Sections 5 and 6 the equations for the transvers and longi- 
tudinal probability density functions are derived and solved. Summary and discussion are 
presented in Section 7. 

Now we would like to recall some well-known properties of velocity correlation functions 
in incompressible fluids, needed below. Consider two points x and x' and define r = x — x'. 
Assuming that the x-axis is paralel to the displacement vector r, one can find that in the 
two-dimensional flow d = 2 for the separation r in the inertial range [5]- [7]: 

! d r r d+1 S 3 = ^£ (9) 



r d + l ' d 

giving 
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S 3 = (A«)3 = (u(x') - u(x)r « (10) 



and 



Si = (Av) 3 = (v(x') - v(x)) 3 ^0 (11) 



where u and v are the components of velocity field paralel and perpendicular to the x- 
axis (vector r). The relations (9)-(ll) resulting from equations of motion (3) are dynamic 
properties of the velocity field. Kinematics also gives something interesting: 



1 d. d -i 



r d 1 dr 

and in two dimensions we have: 



— TJ d - 1 S2 = (d-l)S^(Avy (12) 



1 d 



S 3t = Au(Avy = --S 3 (13) 
3 dr 
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2 Equation for Generating Function 

We consider the iV-point generating function: 

Z =< e Vv(xi) > (14) 

where the vectors x ; define the positions of the points denoted 1 < % < N. Using the 
incompressibility condition, the equation for Z can be written: 

with 

If = E < Aj • f (xj)e A ^^) > (16) 
j 

Ip = -Y J \<e XM ^^> (17) 

The dissipation contributions have been neglected here as irrelevant. 

In what follows we will be mainly interested in the probability density function of the 
two-point velocity differences which is ontained from (7)-(10), setting Ai + A 2 = (see Ref. 
[8] and the theory developed below), so that 

Z =< exp(\ ■ U) > (18) 

where 



U = u(x') - u(x) = Au (19) 

The moments of the two-point velocity differences which in homogeneous and isotropic turbu- 
lence can depend only on the absolute values of two vectors (velocity difference v(x') — v(x) 
and displacement r = x' — x) and the angle 9 between them with 6 = n/2 and 9 = corre- 
sponding to transverse and longitudinal structure functions, respectively. It is easy to show 
[5]- [6] that the general form of the second-order structure function in the inertial range is: 

S2(r,9) = 2 -^D LL (r)(l - ^™s 2 {9)) (20) 
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with Dii(r) =< (u(x) — u(x + r)) 2 >. More involved relation can be written for the fourth- 
order moment: 

S±{r, 9) = D LLLL {r)cos\Q) - 3D LLNN (r)sin 2 (29) + D NNNN (r)sin 2 (9) (21) 

where Dllnn =< (v(x) — v(x + r)) 2 (u(x) — u(x + r)) 2 > and v and u are the components 
of the velocity field perpendicular and parallel to the rr-axis, respectively. In general, in the 
llimit cos{6) = s — > ±1, corresponding to the moments of the longitudinal velocity differences 
S n (r, s) — > S n (r)cos n (9). This means that in this limit Z(X,r, s) — > Z(As,r) = Z(X x ,r). The 
generating function can depend only on three variables: 

A • r / 

Vi = r 5 ^2 = — ^ = Acos(6»); ?7 3 = \j X 2 - i] 2 ; 

In these variables: 

Zt + feA + —d V2 + + {2 ~ d)v2 d m - ^d 2 m ]Z = I f + I p (22) 

f r rr/z r 

where 

I P = Xi< (d 2 , iP (2) - «9 lilP (l)) e A - u > (23) 

and 

I f = (v 2 2+vl)P(l-Cos(k f r))Z (24) 
where, to simplify notation we set d i>a = -^-^ and v(i) = u(xi). 

In two dimensions the equation for the generating function becomes with P = 1 (the subscript 
o is omitted hereafter): 

[d m d m + -«9, 2 + - - (% 2 + % 2 )]^ = / P (25) 

r r o V2 or) 3 r or] 3 2 

The generating function can be written as: 

Z =< e mAu+mAv > (26) 

so that any correlation function 
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< (Au) n Av) m >= ^^Z( V2 = m = 0) (27) 
Neglecting the pressure term I p and differentiating (25) once over i] 2 we obtain immediately 

1 d 

r dr 

Second differentiation (again neglecting I p ) gives: 



-— rS 2 = Si (28) 
r dr 



-^-rS 3 --S 3t -2 = (29) 
r dr r 

Combined with (13) this expression gives 

\±r 3 S 3 -6 = (30) 
r 6 dr 

which is nothing but the Kolmogorov relation, derived in 2d without contributions from 
the pressure terms. It follows from (25) that it is reasonable to look for a scaling solution 

Z(m,V3,r) = Z(X 2 ,X 3 ) where X { = ^A. 

3 Polyakov's theory of Burgers turbulence 

The dissipation-generated contributions 0(vV 2 uiU~) ^ in the limit v — > 0. This is a 



consequence of the ultra-violet singularity V 2 Ui(x)uj(x + r) — > oo when r — > making 
the theory (the closure problem) extremely difficult. The expression for this "dissipation 
anomaly", part of the equation for the generating function, was developed by Polyakov 
for the problem of the one-dimensional Burgers equation stirred by the random force [8]. 
Theory of two-dimensional turbulence is free from the troubles coming from the ultra-violet 
(dissipation) singularities. Still, here we review some of the aspects of Polyakov's theory 
which we believe are of general interest and which will be most helpful below. Polyakov 
considered a one-dimensional problem [8]: 

u t + uu x = f + vu xx (31) 
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where the random force is defined by the correlation function 



f(x,t)f(x + r,t') = K (r)6(t-t) (32) 
The equation for generating function, analogous to (14), is written readily: 

where 

D = u\j < u''^-, t)e Afcu(xfc>t) > (34) 

In the limit — > the force correlation function K{rij) = 0(1 — rf-) which imposes scaling 
properties of the velocity correlation functions. In general, the generating function depends 
on both velocity differences U- = Au = u(x{) — u(xj) and sums U+ = u(xi) + uiXj) which 
makes the problem very difficult. Defining Galilean invariance as independence of the cor- 
relation functions on "non-universal" single-point ul ms = u 2 , Polyakov assumed that if all 
\U-\ « u rms then U- and U+ are statistically independent and ^2Xi — 0. In this case (see 
(8)), introducing ji = \ 2 — Xi and the two-point generating function 

Z(n) =< e^ u > (35) 
the equation for Z reads in a steady state: 

where 

D = iiv < (u"(x + r) - u"{x))e^ u > (37) 

It is clear that the 0(r 2 ) forcing term imposes the scaling variable £ = /ir and Z = F(fj,r) 
where F is a solution of the following equation: 

£F" -F' + i 2 F = D (38) 
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The problem is in evaluation of the dissipation contribution D. 

On the first glance one can attempt to neglect D and solve the resulting equation. This is 
not so simple, however. The Laplace transform of (38) gives an equation for the probability 
density P = = ±$(X) 

$" + x 2 $' + x$ = 

Introducing 

X 3 

$ = Exp( — — )* (39) 

gives 

X 2 

*" = (—- 2X)^ (40) 

which is the Schrodinger equation for a particle in a potential £/(X) = X 4 /4-2X not having 
any positive solutions. 

The positivity of the probability density is a severe constraint on a possible solution of the 
equation of motion. That is where the dissipation contribution D comes to the rescue. 
Polyakov proposed a self-consistent conjecture about the structure of the dissipation term 

D=(- + a)Z (41) 

modifying the potential in the Schrodinger equation with the coefficients b and a chosen to 
produce the zero-energy ground state corresponding to positive PDF. According to Ref.(8) 
this expression is the only one satisfying the galileo invariance of the small-scale dynamics. 
The fact that the one or multi-dimensional advection contributions to the equation for the 
generating function do not lead to positive solutions for the PDF is a general phenomenon 
(see below). The importance of Polyakov's theory is, among other things, in realization that 
the dynamic closures for the remaining terms must remove this problem. This dramatically 
narrows the allowed classes of closures. Thus, the expressions for, D or the pressure terms 
(see below), combined with advective contributions to equation for Z can be correct only and 
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only if they lead to positive solutions for the PDF's in the entire range where \Au\ << u rms 
and r << Lj. 

4 Physical Considerations 

The problem of two-dimensional turbulence is simlified by the fact that the dissipation 
contributions are irrelevant on the scales I » If we are interested in. Moreover, since 
u rms grows with time, the statistically steady small-scale velocity differences U- = Au with 
r « L(t) must be decoupled from U + in (25). This means that the terms 

(Au) n (Av) m (42) 
can eneter the equation for P{ Au, r) while the ones, involving 

(Au) n (Av) m Ul (43) 

cannot. In principle, it can happen that the f/_£/ + -correlation functions can sum up into 
something time-independent. However, at present we discard this bizarre possibility. 
Next, the pressure gradients 

Vp(x + r) - Vp(x) (44) 

appearing in the equation (22)-(24) for Z involve integrals over entire space. It is clear that, 
if the steady state exists, the large- scale contribution to the pressure integrals, depending 
on L = L{t) cannot contribute to the small-scale steady-state dynamics, described by (25). 
That is why the pressure contributions to I p (23) must depend exclusively on the local scale 
r. This leads us to an assumption that the pressure gradients in (23) are local in a sense the 
they can be expressed in terms of the velocity field at the points x and x + r. The application 
of these considerations are presented below. 

The theory of Burgers turbulence, dealt with the "universal" part of the dynamics, i.e. 
with the moments of velocity difference S n with n < 1. The theory of two-dimensional 
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turbulence, we are interested in, must produce the moments with n < oo and that is why 
the algebtaic expressions for the PDF's, characteristic of Burgers dynamics, are irrelevant. 
In addition, we expect the small-scale dynamics in 2d to be independent on the forcing 
function. This makes this problem very different. 



5 Transverse Structure Functions 

Unlike the probability density function for the longitudinal velocity differences P(Au, r), the 
transverse velocity difference probability density is symmetric, i.e. P(Av,r) = P(—Av,r). 
We are interested in the equation (25) in the limit r] 2 — > 0. Let us first discuss some of the 
general properties of incompressible turbulence. Consider the forcing fucntion 



f(x,y) = (f x (x,y),0) 
In this case the equation (25) is: 



l a 'A + ; a - + 7a^-7^-* 2 ' z - / ' < 45 > 

Then, setting 772 = removes all information about forcing function from the equation of 
motion. Based on our general intuition and numerical data we know that two flows strirred 
by a one-component or by a two-component (statistically isotropic) forcing function are 
identical at the scales / >> //, provided the total fluxes generated by these forcing functions 
are equal. This happens due to pressure terms 



Ap = -ViVjViVj 

effectively mixing various components of the velocity field. This universality, i.e. indepen- 
dence of the small-scale turbulence on the symmetries of the forcing, enables us to write an 
expression for the I p contribution to (25). 

According to considerations, presented in a previous section, the pressure gradients in the 
equation (25) are local and their dynamic role is in mixing various components of velocity 
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field. Thus the only contribution to I p , not vanishing in the limit 772 — > 0, can be estimated 
as: 

b— < AuAve V2Au+r > 3Av >=b——< Ave mAu+mAv > ^ 
r r OT]2 

Using a theorem (see Frisch (8), for example) that for the random gaussian process £ (see 
below) 

< mo >=e< ^P- > (47) 

we derive in the limit 77 2 ^ 



2 (At0 2 az 3 



''"W-T-ei (48) 

Substituting this into (25) and integrating over 77 2 gives in the limit 77 2 — > 0: 

dZ 3 Z 3 i] 3 dZ 3 7 2 

-o- + — + —75 + to(m) = T (V3) 49 

or r r or] 3 r z 

where 7 is undetermined parameter and an arbitrary function 

T(r ]3 ) = Z 3 /r + Q(r ]3 ) 

with 

-Q(rj 3 ) = lim V2 ^ ril J Z(rj 2: rj 3: r)drj 2 

is chosen to satisfy a trivial constraint -£3(773 = 0,r) = 1 and the above mentioned universal- 
ity. 

This gives: 

dZ 3 m dZ 3 7 , , 

1T + -o- i --r^ = 50 

where Z 3 = Z(i] 2 = 0,773). This equation is invariant under 773 — > —773 - transformation. It is 
important that the 0(rjf) contribution to (50) comes from the pressure term but not from 
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the forcing, present in the original equation (25). Seeking a solution to this equation in a 
scaling form Z 3 (rj S: r) = Z(r] 3 r^) = Z(X) gives: 

AX 

—Z x = 7 X 2 Z (51) 

and 

Z = Exp^lrl) (52) 

This generating function corresponds to the gaussian distribution of transverse velocity 
differences P(Av) with the second-order structure function 

3 7. 



S*(r) = (Avf = ^r* (53) 

The equation (50) correseponds to a one-dimensional linear Langevin equation for "ve- 
locity field" V = v/{Pr^\ 

v T (x) = —v(x) + (f>(x, t) (54) 

2 1 

where r oc tr~*Pz and the non-local gaussian "universal" forcing (f>{x, r), generated by the 
nonlinearity of the original equation is defined by the correlation function 



(f)(k, r)(f)(k', t') oc 5{k + k')5{r - r') (55) 
The generating function for the field V is 

z =< e xv > 

2 1 

Since r oc tr~* and V oc vr 5 this equation is strongly non-local. It becomes local, however 
in the wave-number space. This will be discussed later. 

Now we can attempt to justify the relation (46). According to (23) and taking into account 
that the rr-axis is paralel to the displacement r in the limit f]2 — > 

I P ~ V3 < (9 y p(0) - d y 'p(r))Exp(r] 3 Av + r] 2 Au) > 
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where 

d y p(0)-d y ,p(r) = J ky(l-e^n[^u(q)u(k-q) + iv(q)v(k-q) + , ^u(q)v(k~q)}d 2 kd 2 q 
and the exponent is expressed simply as: 

eV3 Av+ V2 Au = e iQ^ v ^ d 2 Q +rj2 j( 1 _ e ^) u {^Q)d 2 Q} 

It will be come clear below that transverse velocity differences Av obey gaussian statistics and 
the longitudinal ones Am are very close to gaussian. Then, substituting the above expressions 
into I p and expanding the exponent we generate an infinite series involving various products 
of w(g)'s and t>(g)'s. In case of an incompressible, statistically isotropic gaussian velocity 
field, we are dealing with, these products are split into pairs: 

< VibMQ) >oc q-lfa - ^)5(q + Q) 

The k y integration is carried over the interval — oo < k < oo and in the isotropic case we are 
dealing with the only non-zero terms are those involving even powers of k y . These terms are 
generated by the expansion of 

They however, , being 0(i] 2 ), disappear in the limit r] 2 — > 0. Thus: 

I P = mj d 2 kd 2 qky(l-e ik ^)^ < M (g)t;(A;-g) J Ba;j9(r ? 3y'(l-e^ r )^(Q)^Q+^/(l-e^ r )«(Q)^Q) > 
where the 0{rj 2 ) contribution to the exponent is temporarily kept to make the transformation 

de mAu 
Aue r ' 2 = % 

07]2 

to (46) possible. Only after that we set r] 2 = 0. This proves that the only contribution to the 
equation for the probability density function comes from the O(AuAv) mixing components, 
involved in the pressure gradients. This relation justifies the estimate (46). 
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6 Longitudinal Velocity Differences 

The remarkable fact that in the limit i]2 — > all contributions to the equation (25) contain 
t^- enables separation of variables: integrating the resulting equation over 772 gives the closed 
equation for Z 3 (r] 3 ). The corresponding dynamic equation is linear, meaning that transverse 
velocity fluctuations do not directly contribute to the energy transfer between different scales. 
This effect is possible only in 2d where the 0((d — 2)^-) enstrophy production term in (22), 
not containing is equal to zero. This simplification, combined with locality of the 
pressure-gradient effects, allowed us to derive a closed- form expression for Z3. 

The role of pressure in the dynamics of tranverse components of velocity field is mainly 
restricted to control of the "energy redistribution" neccessary for generation of isotropic 
and incompressible velocity field. The longitudinal field dynamics are much more involved. 
The advection (pressure excluding) part of non-linearity tends to produce large gradients of 
velocity field ("shock generation" using the Burgers equation phenomenology), manifesting 
itself in creation of a constant energy flux in the wave-number space. Pressure is the only 
factor preventing the shock formation. 

Interested in the longitudinal correlation functions we set r] 3 = 0. Then, the term in (25) 

f = f < (Avfe^ *&Z 2 + 0{rtl; % 2 ; rjfa) (56) 

The last relation is accurate since substituting this into (25), differentiating once over r\2 and 
setting both r] 3 = r\2 = gives: 

18 s -_4 = fl^ ,57) 
r or r 3 orj 2 

Since S2 (r) = A 2 ra this equation gives: 

3 01] 2 

which, according to (12) is exact since 9J g^'°- > = (see below). 

Let us consider some general properties of the pressure term I p in the limit 773 — ^ 0. We 
have: 
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dv(2) dv(l) 

I P ~ri2< (-^ - -^)Exp(i l2 Au + r] 3 Av) > (59) 

OX2 OX\ 

Expanding the exponent and recalling that the isotropic and incompressible turbulence Am = 



Av = and p(x)vi(x') = 0, we conclude that 

h « m < - ^-^Au + m Av) 2 + ... >= 0{arf 2 + /377I773 + (60) 

It is clear that the relation (48), derived above for the case of gaussian statistics, satisfied 
this general property of the flow. Thus when 773 — * 0, we approximate 



I p w crf 2 Z + G (61) 

where c is a yet undetermined constant and G denotes the contributions to I p , properly 
modifying numerical coefficients in the equation (25). The presence of the 0(r/|) distinguishes 
this equation from the one for transverse PDF considered in the previous section. There the 
assumed role of pressure was limited to the mixing of various components of velocity field. 
That is why all we accounted for was O(AvAu) contributions to pressure. Here, in addition 
we also consider 0(7)2) contributions, responsible for prevention of the shock formation. The 
resulting equation is: 

r 3 Z 2 - ^tA\^Z 2 - 3r] 2 2 Z 2 - cq\Z 2 = (62) 



r 3 drj 2 dr 5 r 3 

The Laplace transform of gives equation for the probability density P(Au, r): 

1 d 11 A 1 

cPuuu - 3P UV + -^UP + —=-^Pu = (63) 

Seeking solution in a scaling form (the parameter c will be determined below) 

P(U,r) = \F&) (64) 
r 3 r 3 

we obtain 

cF xxx - 3F XX + (b- j)F x + ^xF = (65) 
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Where b = ^A 2 . All, but one, term in (65) change sign wnen x — > —x. The 0(F XX ) 
symmetry-breaking contribution is neccessary for existence of the non-zero energy flux. As- 
suming for a time being that, in accord with numerical and physical experiments, that the 
flux is small (see relation (2)), we first neglect the 0(F XX )- contribution, find solution and 
then take it into account perturbatively. The equation is: 

cF: xx + (b-^-)F: + ^xF° = (66) 

with solution: 

F° = (67) 

A 2 

where c = If A 2 » 1, then the neglected F xx = 0(1/A 2 ) term is small. This means 
that the odd-order moments, computed with the PDF, which is a solution of (65), must be 
small in a sense defined by the relation (2). At the same time the even-order moments must 
be close to the gaussian ones. 

Analytic solution of (65) is difficult. However, one can evaluate all moments % = A n in 
terms of only one parameter A 2 : 

3 A 2 11 

S n+i = — T7t( — ^ n ( n _ l )( n ~ 2 ) S n-3 ~ 3n(n - l)S n - 2 ~ —A 2 nS n - 1 ) (68) 

n + ID 6 6 

This relation gives: A 1 = 0; A 3 = 3/2; A 4 = 3; A 5 = 12A3A 2 ; A 6 = 15Al - 3Q;A 7 = 
37.71A 4 etc. These numbers can be tested in numerical experiments. The one-loop renor- 
malized perturbation expansions give A 2 ps 10, while numerical simulations are consistent 
with A 2 ~ 12. Keeping these numbers in mind, it follows from (68) that the accurate mea- 
surements of the odd-order moments is the only way to verify predictions of the present 
theory. The deviations of the even-order moments from the gaussian ones are too small to 
be detected by both physical and numerical experiments. It can be checked that the ratios 

S 2n +i 

$211+1 2n + l 

C 2n 

D 2n 
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vary in the interval 0.04 — 0.1 for 2 < n < 10 and A2 ~ 10. With A2 ~ 12 these numbers 
decrease even more. 

7 Summary and Conclusions 

The experimentally observed gaussian or very close to it statistics of transverse velocity 
differences was extremely puzzling since, on the first glance, it is incompatible with the 
non-trivial Kolmogorov scaling, resulting from strong non-linearity of the problem. The 
most surpring and interesting result, derived in this paper, is that due to the symmetries 
of the problem the equation, governing probability density function of transverse velocity 
differences, has one derivative less than the one corresponding to the longitudinal differences. 
This means, in turn, that transverse components of velocity field are governed by a non- 
local linear, equation, driven by a universal, non-local, solution-depending gaussian force. 
This reduction, resembling the super- symmetry effects in field theory, is surprising if not 
miraculous. The non-local equation in the physical space, obtained above, corresponds to 
the Langevin equation in the Fourier space: 

v t (k) + c v Phh = f R (k,t) (69) 
where c u is an amplitude of "effective" (turbulent) viscosity and 

f R (k,t)f R (k',t') oc k~ 1 S(k + k')5(t - (70) 

used in [9]- [10] in the renormalization group treatments of fluid turbulence. 
The irrelevance of the dissipation terms in two-dimensional turbulence makes the problem 
much more tractable than its three-dimensional counterpart. Still, in order to close the equa- 
tions for probability density of velocity field one needs an expression for the pressure con- 
tributions. The situation is even more simplified by the fact that the large-scale-dominated 
single-point variables are time-dependent and must decouple from the steady-state small- 
scale dynamics. That is why one can use an assumption about locality of the pressure 
gradient effects leaving only the mixing O(AuAv) contributions to the two-point pressure 
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difference. It can be tested by a mere accounting that all other contributions to the expres- 
sion for I p involve one or more CT+'s and leading to the time-dependent result. This means 
that they must disappear from the steady state equations (25) and (45). The range of pos- 
sible models for pressure is narrowed by a few dynamic and kinematic constraints and by 
the fact that the resulting equation must give positive solution. A simple calculation shows 
that the model for the pressure gradient terms, introduced in this paper, is consistent with 
the derived gaussian statistics. 

The equations for PDF of longitudinal velocity differences do not correspond to linear 
dynamics. Still, the derived solution only slightly deviates from gaussian. This is possible 
due to the relative smallness of the energy flux in two dimensions. 

The results presented here seem to agree with both physical and numerical experiments. 
The obtained close-to-gaussian statistics justifies various one- loop renormalized perturbation 

expansions ^1V1I1^ -^2 ~ 

10 — 12. Using this number we realize that it is extremely difficult to 
experimentally detect deviations from the gaussian statistics. Still, some fine details of the 
present theory, related to the pressure gradient-velocity correlation functions can be tested 
numerically. In addition, measurements of a few odd-order moments can shed some light on 
validity of the present theory. 

The equations and solution presented here leave one question unanswered: are these the 
solutions or not? Our experience with the Burgers and 2d Navier-Stokes equations teach us 
that it is very difficult to find a self-consistent closure leading to the positive solution for the 
PDF's. Stretching this statement a bit we feel that a closure, satisfying dynamic constraints 
and leading to a a plausable solution has a great chances to be correct. 
Absense of intermittency in a steady-state developing inertial range discovered in two- 
dimensional turbulence [2]- [4] seems to be a general phenomenon observed in a drift-wave 
turbulence [11] and in a one-dimensional model of a passive scalar advected by a compress- 
ible velocity field [12]. These observations support our understanding of intermittency as a 
phenomenon originating from interaction of the large and small-scale velocity fluctuations. 
In a developing statistically steady inertial range, were the integral scale is strongly time- 
dependent, these interactions must be small for the small-scale steady state to exist. At the 
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later stages the finite size effects, destroying time-independence of the small scale dynamics, 
lead to formation of coherent structures and new dynamic phenomena which are beyond the 
scope of the present theory 
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